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形 和 斥 才 分 级 依赖 人 工 经 验 判 断 的 现状 ， 本 研究 提出 一 种 采用 轮廓 坐标 系 转换 拟 合 、 
果 形 特征 提取 结合 方向 角 补 偿 算法 检测 柚子 纵 、 横 径 尺 寸 并 基于 果 形 指数 对 柚子 形状 缺陷 进行 判断 的 方法 。 


以 CMOS 相 机 、 点 阵 式 LED 光源 、 平 面 镜 、 计 算 机 、 箱 体 和 支架 搭建 图 像 采 集 装 置 ， 获 取 168 个 不 同 尺 寸 与 


形状 等 级 的 沙 昌 


昌 柚 样本 全 表面 图 像 数 据 。 选 择 G-B 分 量 灰 度 图 像 进行 去 品 与 分 割 ， 利 用 Laplacian 算 子 边 缘 


检测 算法 提取 果实 的 边缘 像素 ,采用 多 项 式 拟 合 方式 完成 直角 坐标 向 极 坐 标的 转换 从 而 简化 果 形 描述 ， 利 
用 特征 点 极 角 值 补偿 样本 纵横 径 的 随机 方向 ， 继 而 区 别 类 球形 和 类 梨 形 两 种 类 型 计算 柚子 的 纵 径 和 横 径 。 


以 广东 梅州 沙田 


柚 为 对 象 进行 试验 ， 结 果 表 明 ， 利 用 轮廓 坐标 系 转换 拟 合 、 果 形 特征 提取 结合 方向 角 补 偿 


算法 的 方法 检测 柚子 纵 径 的 平均 绝对 误差 、 最 大 绝对 误差 和 平均 相对 误差 分 别 为 2.23 mm、7.39 mm 和 1.6%， 


横 径 的 平均 绝对 


误差 、 最 大 绝对 误差 和 平均 相对 误差 分 别 为 2.21 mm、7.66 mm 和 1.4%。 从 柚子 轮廓 极 坐标 


的 拟 合 函数 中 提取 3 个 峰值 高 度 、3 个 波峰 宽度 和 1 个 波 谷 值 差 值 7 个 特征 值 ， 利 用 BP 神 经 网 络 算法 建立 柚 


果 形 判别 模型 并 用 独立 验证 集 进行 验证 ， 形 状 判 别 的 总 识别 率 为 83.7%。 本 方法 能 为 柚子 尺寸 和 形状 的 自 


动 化 检测 与 分 级 提供 快速 无 损 方法 。 
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1 引言 


柚子 (Citrus maxima (Burm) Merr.) 是 东 
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态势 "。 作 为 一 种 体积 较 大 的 内 部 非 均 质 水 果 ， 
柚子 果实 形状 因 品 种 和 生长 条 件 不 同 呈 梨 形 、 近 
梨 形 、 球 形 等 ， 其 果 形 特征 通常 用 果 带 、 柚 颁 、 
果 户 、 杰 道 和 果 顶 等 部 位 的 尺寸 和 部 位 之 间 比 例 
关系 进行 描述 ， 如 图 1 所 示 。 
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图 1 柚子 果 形 特征 


Fig. 1 Shape features ofpomelo fruit 


常见 的 柚子 果 形 缺 隐 有 柚 颈 过 长 、 不 对 称 、 
球形 和 不 符合 品种 特征 等 。 因 此 ， 仅 以 重量 值 或 
最 大 尺寸 值 为 指标 进行 柚子 分 级 均 不 能 满足 按 质 
论 价 的 要 求 。 当 前 生产 中 部 分 加 工 企业 借助 通用 
型 果品 分 级 设备 依据 重量 或 过 租 原 理 任 有 果实 大 小 
区 分 档次 ， 不 能 完全 符合 柚子 分 级 的 特性 要 求 ; 
果园 产地 分 级 主要 依赖 人 工 ， 依 据 柚子 大 小 和 形 
状 凭 经 验 粗 略 处 理 ， 不 仅 劳 动量 大 、 效 率 低下 ， 
而 且 判 断 标 准 无 法 统一 。 加 之 中 国 国内 柚子 多 在 
山地 丘陵 地 区 的 中 小 果园 种 植 ， 劳 动力 短缺 日 益 
明显 ， 产 后 综合 分 级 技术 缺乏 等 ， 这 些 因素 严重 
制约 了 柚子 产业 的 效益 提升 和 规模 发 展 。 因 此 ， 
研究 一 种 能 够 根据 柚子 果 形 特点 进行 自动 化 检测 
分 级 的 方法 ， 对 于 柚子 生产 具有 重要 意义 ”。 

机 器 视觉 技术 在 获取 待 测 对 象 外 部 品质 参数 
中 具有 独特 优势 ， 能 够 快速 无 损 地 获取 果蔬 对 象 
的 尺寸 、 颜 色 和 表面 缺陷 等 参数 “， 结 合 特征 信 
息 提取 和 人 处 理 算法 鉴别 果蔬 品质 。 例 如 ， 王 红军 
等 ”探讨 了 利用 机 器 视觉 结合 主 成 分 分 析 法 和 
多 元 线性 回归 方法 ， 筛 选 出 权重 系数 大 的 图 像 特 
征 数 据 ， 建 立 了 马铃薯 质 量 和 形状 分 级 预测 模 
型 。 王 风云 等 ”选用 双 苑 蘑 菇 菌 盖 直径 作为 尺 
才 分 级 的 特征 参数 ， 结 合 全 局 阔 值 分 割 法 与 寻找 
Te BUA. 、Canny 算 子 、 闭 运算 、 分 水 
岭 算 法 去 除 图 像 中 阴影 部 分 和 柄 部 干扰 ， 采 用 最 
小 外 接 矩 形 法 求 得 了 菌 盖 直径 。 黄 展 和 费 继 友 
提出 了 改进 的 三 层 Canny 边缘 检测 算法 来 提取 苹 
果 轮 廓 ， 并 根据 标准 矩形 工件 斥 十 来 标定 像素 当 


量 ， 进 而 测量 果 径 、 面 积 等 特征 。Roscher 等 " 
通过 采集 自然 环境 下 的 葡萄 图 像 ， 采 用 Hough 检 
测 圆 的 方法 检测 所 有 可 能 的 类 球 型 葡萄 ， 通 过 颜 
色 、 纹 理 等 信息 设计 分 类 避 ， 检 测 出 葡萄 的 数量 
与 果 径 大 小 ， 与 人 工 测量 相 比 ， 测 得 的 果 径 平均 
差 值 约 为 1 mm, 平均 浆果 大 小 与 人 工 测量 之 间 
的 相关 性 为 0.88， 同 时 可 有 效 降低 检测 错误 率 。 
Oo 等 ”将 草 每 形状 简化 为 风筝 (kite)， 从 水 果 
的 边界 中 识别 出 上 、 下 、 左 、 右 四 个 顶点 及 任意 
两 个 癌 量 之 间 的 夹 角 ， 以 果实 项 点 与 顶点 之 间 的 
垂直 距离 作为 草 每 的 长 度 ， 以 两 侧 点 之 间 的 距离 
作为 果实 的 直径 。Arendse 等 ”利用 X 射 线 计算 
机 汤 层 扫描 技术 将 石榴 果实 生成 二 维 放 射 图 像 ， 
通过 三 维 重建 ， 实 现 对 果实 长 度 、 直 径 、 果 皮 厚 
度 等 参数 的 估计 。 上 述 研究 都 是 利用 果实 直径 、 
长 度 、 面 积 等 参数 对 果实 形状 进行 估计 ， 而 果实 
的 果 形 特征 一 般 凭借 人 工 的 经 验 去 判别 ， 果 实 类 
型 不 同 判 别 方法 也 有 很 大 差别 ， 故 而 检测 果 形 的 
方法 实现 难度 较 大 ， 相 关 人 研究 也 相对 较 少 。 目 前 
较 篆 见 的 果 形 特征 判别 方法 有 提取 轮 廊 线 和 三 维 
重建 法 等 。 罗 陆 锋 等 …"' 利用 改进 K-means 聚 类 
方法 对 葡萄 图 像 进行 分 割 ， 提 取 轮 廓 边缘 轮廓 和 
左右 轮廓 的 类 圆 中 心 ， 以 中 心 点 为 原点 建立 基于 
轮廓 分 析 的 二 贴 葡萄 串 分 界线 几何 求解 和 计算 模 
型 ， 通 过 求解 两 玛 由 葡萄 边缘 轮廓 交界 处 拐点 来 
实现 对 双 串 三 贴 和 葡萄 目标 分 别提 取 。 有 吴 明 清 
等 1) 通过 提取 红 束 的 二 维 轮廓 特征 ， 由 二 维 图 
像 的 轮廓 构建 三 维 多 轮 廊 模 型， 以 此 检测 红枣 的 
形状 。Olatunji 等 "” 利用 条 件 生 成 对 抗 网 络 对 生 
成 的 大 量 猕猴 桃 形 状 信息 进行 训练 ， 来 真实 地 重 
建 封闭 和 完整 的 猕猴 桃 表面 。 

鉴于 柚子 不 完全 对 称 的 特性 ， 并 且 现 有 的 沙 
田 查分 级 标准 对 果 形 的 描述 过 于 笼统 ， 缺 乏 针对 
具体 果 形 参数 的 分 类 标准 ， 本 研究 提出 通过 获取 
柚子 果实 全 表面 图 像 信 息 ， 采 用 轮廓 坐标 系 转换 
拟 合 、 果 形 特征 提取 结合 方向 角 补 偿 算法 检测 机 
子 纵横 径 尺 寸 ， 并 基于 果 形 指数 对 柚子 形状 缺陷 
进行 判断 的 方法 ， 以 期 为 柚子 果 形 检测 与 分 级 提 
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供 一 种 快速 无 损 的 解决 方案 。 
2 试验 材料 与 装置 


2.1 试验 样本 获取 


在 广东 省 梅州 农业 实验 站 和 周边 果园 选 购 采 
收 期 沙田 柚 类 球形 样本 81 个、 类 梨 形 样本 54 个 
AUER GR CHR VD FAA RIE TIE (+) 样本 
33 个 ， 共 计 168 个 试验 图 像 作 为 采集 用 果实 样 
本 。 对 每 个 样本 进行 尺寸 测量 与 形状 判定 后 进行 
编号 ， 并 进行 图 像 数据 采集 。 借 助 可 在 水 平 导轨 
上 滑动 的 两 块 直角 夹板 将 果实 按 待 测 纵 轴 或 模 轴 
方向 夹 紧 后 ， 用 数 显 式 游标 卡尺 对 每 个 样本 测量 
2 次 取 平 均值 作为 该 样本 的 纵 径 或 横 径 值 。 样 本 
形状 判定 依据 沙田 柚 行业 标准 NY/T 868- 
2004'”， 通 过 人 工 观 察 与 测量 沙田 柚 样本 柚 颂 
长 度 、 纵 径 与 横 径 比例 和 形状 规则 度 ， 将 试验 样 
本 分 为 3 类 ,样本 外 观 如 图 2 所 示 。 


(DERE: (b) 第 一 类 缺陷 : (c) 第 二 类 缺陷 : 
标准 果 形 柚 颈 或 果 肩 过 长 类 球形 


图 2 3 类 柚子 样本 原始 图 像 


Fig. 2 Original images ofthree pomelo types 
2.2 图 像 采 集 系 统 


利用 机 器 视觉 和 图 像 处 理 技术 对 样本 图 像 的 
果 形 特征 进行 提取 和 分 析 具 有 非 破坏 性 、 互 容 性 
高 等 优势 。 图 3 为 自行 搭建 的 图 像 采集 系统 结构 
示意 图 。 该 系统 主要 包括 工业 相机 、 光 源 、 图 像 
采集 箱 体 和 计算 机 等 部 分 。 工 业 相 机 为 CMOS 相 
机 (台湾 显 微 1/2.5 MU500， 中 国 ) ， 最 大 分 辩 率 
2592x1944 像素 ， 镜 头 型 号 XW0612。 图 像 采集 
箱 体 外 部 尺寸 80 cm x60 cmx60 cm， 箱 体内 部 布 
置 防 反光 绒布 ， 消 除 反 射 光 影响 ， 采 用 12 WAY 
LED 贴 片 点 阵 式 圆 形 为 主 光 源 及 10 W (可 调 ) 


贴 片 点 阵 式 辅助 光源 。 箱 体内 设置 相对 放置 的 两 
块 平面 镜 用 于 使 相机 单 次 拍摄 可 获取 类 梨 样 本 的 
全 部 表面 信息 ""。 经 分 析 柚 子 尺寸 范围 和 前 期 
调试 ,平面 镜 尺寸 取 30 cm x 25 cm， 相 邻 端点 间 
距 取 20 cm, 镜面 夹 角 取 120° 以 保证 相机 获取 的 
图 像 能 够 包含 样本 全 部 表面 信息 。 


注 :1. 辅 助 光源 2. 平 面 镜 3.CMOS 相机 4. 主 光源 5. 样 本 放置 位 
图 3 图 像 采集 系统 结构 图 


Fig. 3 Structure diagram of image acquisition system 


图 像 采 集 时 样本 放置 于 样本 位 ， 方 向 随机 。 
利用 Minte Camera Platform 软件 控制 相机 获取 样 
本 图 像 ， 并 经 由 数据 线 传人 计算 机 进行 储存 和 分 
析 。 系 统 所 得 样本 图 像 如 图 4 所 示 。 


图 4 系统 获取 的 原始 样本 图 像 示例 


Fig. 4 A random original image obtained by the system 
3 ”特征 提取 及 形状 检测 


本 人 研究 对 采集 的 柚子 图 像 进行 预 处 理 操 作 之 
后 ， 利 用 极 坐标 系 代替 传统 的 平面 直角 坐标 系 ， 
提取 柚子 轮廓 的 特征 信息 。 将 柚子 按 果 形 分 为 类 
BIE MIRE DAS AHI A tas FU > R UE LAY 
差 值 检测 果实 的 模 径 和 纵 径 ， 通 过 BP 神经 网 络 
建立 柚子 形状 判别 模型 对 果实 进行 分 类 。 图 5 为 
本 研究 技术 路 线 。 
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像 预 处 理 


轮廓 坐标 系 转换 


果实 横 径 、 纵 径 检测 


形状 判别 模型 建立 及 验证 


基于 轮 慷 坐标 系 转换 拟 合 的 
柚子 果 形 检测 分 级 方法 


图 5 特征 提取 及 形状 检测 技术 路 线 
Fig. 5 The technical route of feature extraction 


and shape detection 


3.1 图 像 预 处 理 


由 于 在 拍照 过 程 中 ， 环 境 变化 和 噪声 干扰 均 
会 使 得 拍照 所 得 的 图 像 特征 不 明显 ， 故 需 对 图 像 
进行 滤波 去 噪 、 特 征 增强 等 处 理 手段 ， 实 现 图 像 
特征 突出 的 效果 ,方便 后 续 图 像 识 别 。 比 较 均值 
滤波 、 中 值 滤波 等 图 像 滤波 方法 ， 和 迭代 法 、 最 大 
类 间 访 查 法 等 图 像 分 割 方法 以 及 Canny、Sobel、 
Robert, Laplacian 等 边缘 算 子 提取 算法 对 样本 图 
像 的 处 理 效果 。 综 合 考 虑 轮廓 信息 提取 效果 和 运 
SIE, RAE REBT, ext 
图 像 进 行 高 斯 滤波 平滑 图 像 、 去 除 噪声 ， 在 
RGB 颜色 模型 下 选取 G-B 色差 的 图 像 分 割 特征 
描述 柚子 图 像 "… ， 如 图 6 所 示 。 利 用 最 大 类 间 
TÆ (Otsu) 对 灰 度 网 像 进行 二 值 化 处 理 ， 
利用 中 值 滤波 算法 去 除 二 值 图 干扰 噪声 像素 ， 经 
颜色 国 值 分 割 后 获得 与 背景 分 离 的 目标 果实 图 
像 ， 二 值 图 像 开 运算 保持 了 轮廓 大 小 并 填补 二 值 
图 像 中 的 白色 孔洞 和 细 缝 ， 最 后 借助 Laplacian 
边缘 检测 算 子 提取 目标 图 像 的 有 效 边 缘 信 
息 ""”"， 样本 边缘 提取 结果 图 像 如 图 7 所 示 。 


3.2 轮廓 坐标 系 转换 拟 合 


图 像 预 处 理 后 所 得 样本 轮廓 为 自然 形状 ， 果 
实 纵 轴 方 向 随机 。 由 于 柚子 果 贷 和 果 肩 的 曲线 弧 
度 变 化 和 方向 随机 分 布 ， 使 得 轮廓 曲线 在 平面 直 


图 6 样本 G-B 色差 图 像 


Fig. 6 Image of chromatic aberration G-B 


component of the sample 


图 7 样本 边缘 轮廓 图 像 
Fig. 7 Edge contour image of the sample 

角 坐 标 系 中 的 函数 复杂 且 难 于 确定 纵横 轴 的 方 
向 ， 纵 径 与 横 径 长 度 的 计算 量 大 ， 检 测 与 分 级 的 
误差 增加 。 极 坐标 系 属 二 维 坐标 系统 ， 可 与 平面 
直角 坐标 系 进 行 相 互 转换 。 相 较 于 平面 直角 坐标 
系 ， 极 坐标 系 可 以 简化 对 于 柚子 轮廓 曲线 的 描 
述 ， 故 而 将 轮廓 曲线 函数 从 直角 坐标 系 向 极 坐标 
系 进 行 转换 拟 合 的。 利用 公式 C1) 和 公式 (2) 
获取 柚子 质心 坐标 ， 选 取 直 角 坐 标 系 义 轴 负 方向 
为 极 轴 方 向 ， 极 坐标 系 中 横 轴 表示 极 角 角度 ， 纵 
轴 表 示 极 径 ， 从 而 实现 将 直角 坐标 系 中 的 轮廓 转 
换 至 极 坐标 系 进行 描述 ， 得 到 样本 轮廓 信息 的 极 
坐标 图 小 ， 如 图 8 所 示 。 通 过 极 坐 标 图 中 极 径 最 
大 峰值 可 确定 柚子 的 纵 径 与 水 平方 向 的 偏 移 角 
度 。 为 避免 极 坐标 描述 数据 计算 量 过 大 增加 数据 
处 理 量 和 计算 时 间 ， 采 用 间隔 固定 步 长 选取 像素 
点 法 缩减 数据 量 。 图 9 比较 了 步 长 m 在 [10, 30] 
区 间 内 以 10 个 步 长 为 间隔 取 值 的 效果 ， 可 见 当 m 
为 20 时 ， 有 效 减少 了 计算 量 的 同时 最 大 程度 保 
留 了 柚子 果 形 轮廓 信息 。 
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d 
y= >/s (2) 
i=0 
其 中 > 2 为 质心 点 X 轴 坐标 ; 区 为 质心 点 Y 
轴 坐 标 ; 5 为 轮廓 总 像素 数 ， 个 。 
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图 8 不 同 坐 标 系 中 的 类 梨 形 柚子 样本 轮廓 


Fig. 8 The edges of a pear-like sample in cartesian 


coordinate and polar coordinate systems 
3.3 极 坐 标 曲线 拟 合 


曲线 拟 合 (Curve Fitting) 是 采用 适当 曲线 


w 
a 
So 
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(c) 步 长 m=30 


图 9 步 长 取 10、20 和 30 的 轮廓 描述 
Fig. 9 Contours description based on step lengths 


at 10,20 and 30 


用 方向 角 确 定 柚子 纵 径 方向 ， 再 通过 计算 果实 轮 
廓 图 像 X 方 向 和 YY 方向 的 最 大 距离 来 确定 类 梨 形 
柚子 的 纵 模 径 尺 寸 大 小 ; 对 于 类 球形 柚子 ， 采 用 
最 小 外 接 和 矩形 法 检测 其 尺寸 。 类 梨 形 柚子 果 形 在 


函数 对 符 测 数据 进行 描述 的 方法 ， 并 通过 拟 合 后 
的 结果 对 变量 数据 进行 分 析 处 理 ， 进 而 获得 数据 
间 的 变化 规律 。 本 研究 采用 公式 (3) 所 表示 的 
多 项 式 拟 合 (Polynomial Fitting) 方法 进行 柚子 
轮廓 信息 的 数据 拟 合 处 理 。 拟 合 结果 以 相关 系数 
R 和 均 方 根 误差 RMSE 作为 评价 指标 进行 评价 ， 
其 中 R 越 大 、RMSE 越 小 则 表明 拟 合 结果 越 好 。 
SS (3) 
HP, x AMA AE et x I TES ABI 
输入 变量 系数 ; kA USK RL 
3.4 纵 径 和 横 径 的 图 像 检测 方法 


为 提高 检测 精度 ， 有 效 减 小 误差 ， 将 样本 区 
分 为 类 梨 形 和 类 球形 两 类 。 对 于 类 和 巢 形 柚 子 ， 利 


直角 坐标 系 和 极 坐 标 系 中 的 轮廓 如 图 8 所 示 ， 类 
球形 柚子 果 形 在 直角 坐标 系 和 极 坐 标 系 中 的 轮廓 
如 图 10 所 示 。 在 极 坐 标 系 中 可 以 观察 到 ， 极 径 
最 大 峰值 和 第 二 个 峰值 之 差 的 差异 最 为 显著 ， 故 
利用 公式 (4) 计算 极 径 峰 值 差 值 比 C,。 随 机 抽 
取 9 个 类 球形 和 9 个 类 梨 形 样本 的 数据 计算 C; 的 
值 并 用 于 对 比分 析 ， 最 终 确 定 利 用 峰值 差 值 比 C, 
判别 柚子 类 球形 或 类 梨 形 的 分 类 阔 值 为 0.34。 
hi= hy 

h, 

其 中 ，C; 为 极 径 峰值 差 值 比 ; hh 为 极 径 最 大 
E, mm; 思 为 极 径 第 二 个 峰值 ，mm。 

对 类 梨 形 样 本 纵横 径 尺寸 进行 检测 。 首 先 确 
定 其 纵 轴 方向 ， 然 后 再 计算 纵横 径 尺 寸 。 图 11 
为 方向 角 在 -90"~90" 范 围 内 的 柚子 样本 二 值 图 


C, = (4) 
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Fig. 10 The edges of a spherical sample in cartesian 
coordinate and polar coordinate systems 

像 ， 结 合 样本 轮廓 的 极 坐标 图 像 进行 分 析 可 知 ， 
样本 轮廓 最 大 极 径 处 即 为 该 样本 的 果 蒂 位 置 ， 故 
而 最 大 极 径 值 处 所 对 应 的 极 角 角度 即 可 用 于 替代 
子 纵 径 的 方向 角 。 

表 1 列 出 了 对 应 样本 的 方向 角 和 极 角 角度 
值 ， 利 用 极 角 角度 奉 代 实 际 纵 径 方 向 角 的 平均 误 


> 


(a)-52.5° (b)2.2° 


(c)30.2° 


(d)85.3° 
图 11 不 同方 位 柚子 图 像 
Fig. 11 Images of pomelo samples with different 


orientations 


差 为 2.55°* ， 表 明 可 利用 极 角 角度 来 进行 柚子 纵 
径 方 向 的 估算 。 以 极 角 角 度 值 进 行 图 像 旋 转 可 得 
果实 纵 轴 垂 直 向 上 的 轮廓 图 像 ， 计 算 该 轮廓 图 像 
在 Y 方 向 和 X 方 向 的 最 大 差 值 即 为 该 类 和 巢 形 样本 
的 纵 径 和 横 径 尺寸 。 

表 1 不 同方 位 柚子 样本 方向 角 与 极 角 间 的 关系 


Table 1 The relationship between direction angel and polar 


angel ofpomelo with different orientations 


图 像 编 号 方向 角 /(°) fC) 误差 /(°) 
a -52.5 -49.7 -2.8 
b 2.2 0 2.2 
c 30.2 27.4 2.8 
d 85.3 82.9 2.4 


类 球形 柚子 由 于 纵横 径 差 值 较 小 ， 以 方向 角 
确定 纵 轴 方向 的 方法 并 不 适用 ， 故 采用 最 小 外 接 
矩形 法 (Minimum Enclosing Rectangle) 2 进行 
果实 尺寸 检测 。 根 据 柚子 质心 描绘 出 柚子 的 外 接 
矩形， 然后 让 柚子 边缘 围绕 质心 进行 顺 时 针 或 道 
时 针 旋 转 ， 每 旋转 一 次 便 计 算 一 次 随 之 改变 的 和 矩 
形 面积 ， 经 过 360? 的 旋转 后 ， 面 积 最 小 的 矩形 便 
是 此 柚子 的 最 小 外 接 和 矩形 ， 其 最 小 外 接 和 矩形 的 长 
和 宽 即 可 作为 样本 的 纵 径 和 横 径 。 


3.5 形状 判别 模型 建立 


通过 分 析 样 本 轮廓 极 坐标 特性 peh 
纵 径 和 横 径 用 于 特征 描述 。 首 先 判别 纵 径 小 于 
152 mm 或 横 径 小 于 142 mm 的 过 小 等 外 果 ， 然 后 
R (H) 和 横 径 (D) 的 比值 为 果 型 的 判别 
指数 S (S=HD)， 若 S 大 于 1.1， 则 说 明 该 样本 
果 颈 长 度 适宜 ， 果 形 合 格 o 
通过 比较 随机 取样 法 (Random Sampling， 
RS)"”、 肯 纳 德 -斯 通 算 法 (Kennard-Stone， 
KS)" 、 基 于 联合 xy 距离 的 样本 集 划 分 法 
(Sample Set Partitioning Based on Joint x-y Dis- 
tance，SPXY) ”三 种 样本 集 划 分 方法 的 判别 效 
有 果 。RS 法 操作 简单 方便 , 但 由 于 过 程 随机 ， 难 
以 保证 建立 模型 的 稳定 性 ; KS 法 实现 了 校正 集 
在 空间 距离 上 的 均匀 分 布 ， 但 计算 量 过 大 ; 
SPXY 法 与 KS 法 类 似 ， 同 样 利 用 样本 的 欧 氏 距 
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离 进 行 校正 集 的 划分 ， 但 同时 将 x 变量 和 y 变量 
考虑 其 中 ， 可 以 覆盖 多 维 空间 ， 增 强 预 测 模型 效 
果 。 综 合 考虑 ， 选 取 SPXY 法 按照 校正 集 样 本 
119 个 、 验 证 集 样本 49 个 划分 样本 集 。 利 用 BP 
神经 网 络 建立 判别 模型 ， 通 过 不 断 更 新 网 络 中 的 
冰 值 和 权重 ， 实现 最 终 输 出 结果 向 理想 值 的 通 
近 。 通 过 比较 隐 含 层 节 点 数 为 5 一 10 个 的 网 络 结 
构 和 tansig Phi AC. logsig 函数 和 purelin 函数 三 种 
不 同 的 节点 传递 函数 对 BP 神经 网 络 模型 建 模 判 
别 结果 的 影响 ， 确 定 最 优 BP 神 经 网 络 模型 的 建 
立 条 件 的 ， 选 择 隐 含 层 节 点 数 为 6，tansig 函数 
为 节点 传递 函数 建立 判别 模型 。 图 像 处 理 和 模型 
建立 均 利 用 MATLAB R2014b (MathWorks 公司 ， 
美国 ) 实现 。 


4 结果 与 讨论 


4.1 轮廓 坐标 系 转换 拟 合 结果 


利用 公式 (1) 和 (2) 对 样本 轮廓 图 像 由 直 
角 坐 标 系 向 极 坐标 系 转 换 后 ， 根 据 公 式 (3) 进行 
多 项 式 拟 合 。 表 2 列 出 了 取 不 同 多 项 式 最 高 次 数 
时 对 样本 轮廓 曲线 的 拟 合 结果 ， 可 见 上 取 15 时 ， 
相关 系数 R 为 0.9714。k 增 加 到 17 时 ，R 增 加 至 
0.9736， 此 后 上 增加 而 不 变 。 故 取 k 为 17 时 ， 
此 时 RMSE 为 0.0190， 运 算 时 间 需 0.000103 s, 
能 够 反映 样本 轮廓 特征 并 满足 处 理 速度 要 求 。 
表 2 不 同 多 项 式 最 高 次 数 对 样本 轮廓 曲线 的 拟 合 结果 
Table 2 Fitting results of the highest degree of different 


polynomials to the sample contour curve 


多 项 式 拟 合 R RMSE 拟 合 时 间 /s 
k=15 0.9714 0.0198 0.000084 
k=16 0.9713 0.0198 0.000092 
k=17 0.9736 0.0190 0.000103 
k=18 0.9736 0.0190 0.000153 


4.2 柚子 样本 纵 径 和 横 径 的 检测 结果 


利用 数 显 式 游标 卡尺 人 工 对 168 个 柚子 样本 进 
行 尺寸 测量 ， 所 得 纵 径 范围 为 100.10~218.22 mm, 
横 径 范围 为 88.74~161.36 mm。 利 用 1 张 图 像 数 
据 上 的 3 幅 样本 轮廓 进行 图 像 分 析 ， 分 别 计 算 主 


视图 、 左 镜像 、 右 镜像 3 个 方位 柚子 图 像 的 尺寸 
大 小 ， 选 取 3 者 的 均值 作为 该 柚子 样本 最 终 的 尺 
才 大 小 计算 值 ， 检 测 计算 值 与 人 工 测量 值 间 的 各 
项 误差 如 表 3 所 示 。 横 径 的 平均 相对 误差 为 
1.4%， 纵 径 的 平均 相对 误差 为 1.6%。 

RI 样本 纵 、 横 径 检 测 结 果 误 差 


Table 3 Test result errors of the samples vertical and 


horizontal diameter 


最 大 绝对 最 大 相对 ”平均 绝对 ”平均 相对 


检测 指标 
误差 /mm 误差 /% 误差 /mm ” 误差 /% 
横 径 7.66 4.5 2.21 1.4 
纵 径 7.39 4.5 2.23 1.6 


4.3 果 形 特征 模型 建立 及 验证 


对 拟 合 函数 求 二 阶 导 ， 提 取 3 类 果 形 样本 的 
3 个 峰值 高 度 、3 个 波峰 宽度 和 1 个 波 谷 值 差 值 共 
7 个 特征 值 用 以 描述 柚子 轮廓 拟 合 函 数 ， 如 图 12 
所 示 。 


200 
-200 -150 -100 -50 0 50 100 150 200 
极 角 /(”) 


图 12 拟 合 函数 特征 
Fig. 12 Features of fitting function 

样本 总 量 为 168 个 ， 利 用 119 个 校正 集 样本 
建立 BP 神经 网 络 ， 初 始 设 置 网 络 结构 为 7-9-1， 
其 中 包含 7 个 输入 层 节 点 、9 个 隐 含 层 节 点 、1 个 
输出 层 节 点 ， 学 习 速 率 为 0.1， 浆 值 随机 产生 更 
新 。 比 较 分 类 模型 的 判别 正确 率 ，6 个 隐 含 层 
节点 的 网 络 模 型 判别 准确 率 达 84% ， 为 最 高 ; 
tansig PAX. logsig PK 2X Fl purelin 函数 判别 准确 
率 分 别 为 84%、80% 和 52%。 确 定 最 优 的 BP 神 
经 网 络 结构 为 7-6-1， 包含 7 个 输入 层 节 点 、6 个 
隐 含 层 节 点 、1 个 输出 层 节 点 ,最 优 节 点 传递 函 
数 为 tansig 函数 。 隐 合 层 每 一 个 节点 会 对 输入 层 
的 兴奋 有 不 同 的 接收 权重 ， 从 而 更 加 偏向 于 某 种 
识别 模式 。 该 模型 输入 层 与 隐 含 层 的 权重 如 下 ， 
和 矩阵 的 行 表示 隐 含 层 节 点 ， 列 表示 输入 层 节 点 ， 
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| -0.1446 0.4797 -0.7768 0.6710 -0.0766 0.8888  —0.4414| 
0.7279 -1.9122 0.0906 -1.2852 -0.0446 1.9385 1.0107 
1.7425 -1.4097 0.6493 -0.0689 0.8272 1.4207 -0.7080 
-2.1120 -2.3170 1.0909 -0.7888 -0.1764 -1.8285 0.7673 
-0.0392 3.7368 0.6586 -3.1757 0.6753 -0.7993 0.7550 

L 1.4191 -0.3585 -0.7487 0.7783 -1.4139 0.3027 1.0698 | 


经 过 最 优 判 别 ， 隐 含 层 与 输出 层 的 权重 为 
[ -0.6010, 0.6917, 0.3088, -0.5669， -0.8380, 
-2.2320]; ka & Jz bd fa A | -1.9280, -2.2626, 
-0.1841, 0.2974, 1.5112, 1.7438] 7; $f H JE BJE 
为 2.2937， 即 为 柚子 果 形 指数 。 根 据 柚子 实际 果 
形 与 果 形 系数 的 关系 ,定义 柚子 果 形 判别 标准 为 : 
若 输 出 果 形 指数 在 (-%, 1.5) 之 间 ， 则 该 样本 为 
果 形 标准 样本 ; 若 输 出 果 形 指数 在 [1.5, 2.5) 之 
间 ， 则 该 样本 为 柚 颈 或 果 肩 过 长 的 第 一 类 果 形 缺 
陷 样 本 ; 若 输 出 果 形 指数 在 [2.5, +0) 之 间 ， 则 
该 样本 为 缺乏 沙田 柚 果 形 特征 的 第 二 类 类 球形 果 
形 缺 陷 样 本 。 

利用 所 建 模型 对 49 个 验证 集 样本 进行 判别 ， 
其 结果 如 图 13 所 示 。 类 别 1 为 29 个 果 形 标准 的 
柚子 样本 ， 误 检 数 为 1 个 ， 漏 检 数 为 5 个 ; 类 别 
2 为 15 个 果 形 较 规 则 ， 果 形 缺 陷 较 小 ， 柚 颈 或 果 
肩 过 长 的 柚子 样本 ， 误 检 数 为 5 个 ， 漏 检 数 为 2 
个 ; 类 别 3 为 5 个 果 形 缺陷 较 大 ， 类 球形 缺乏 沙 
田 柚 的 果 形 特征 ， 误 检 数 为 2 个 ， 漏 检 数 为 1 
个 。 综 上 ， 计 算得 到 该 模型 总 识别 率 为 83.7%。 
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本 研究 采集 168 个 柚子 样本 的 全 表面 图 像 数 
据 ， 利 用 传统 平面 直角 坐标 系 转换 为 极 坐 标 系 简 
化 柚子 果 形 轮廓 描述 ， 同 时 消除 查 果 方向 位 置 随 
机 放置 导致 的 纵横 径 判别 困难 ， 角 度 校 正 的 平均 
误差 为 2.55°。 试 验 样本 尺寸 检测 的 纵 径 平均 误 
差 为 2.23 mm, 平均 相对 误差 1.6%; 横 径 平均 误 
差 为 2.21 mm, 平均 相对 误差 1.4%。 处 理 时 间 需 
0.000103 s。 

利用 样本 轮廓 的 极 坐 标 拟 合 曲 线 进行 柚 子 果 
形 描述 ， 提 取 7 个 柚子 果 形 特征 建立 BP 神经 网 
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图 13 BP 神经 网 络 预 测 分 类 与 实际 类 别 比 对 图 
Fig.13 Back propagation neural network prediction classifica- 


tion and actual category comparison 


络 判别 模型 ， 确 定 最 佳 的 网 络 结构 为 7-6-1， 节 
点 传递 图 数 为 tansig 国 数 ， 模 型 总 识别 率 为 
83.7%。 研 究 结果 表明 ， 采 用 轮廓 坐标 系 转换 拟 
合 、 果 形 特 征 提取 结合 方向 角 补 偿 算 法 进行 柚子 
纵 、 横 径 尺寸 检测 和 基于 果 形 指数 的 形状 缺陷 判 
断 方法 能 够 实现 基于 机 融 视 觉 的 柚子 果 形 指标 快 
速 无 损 检 测 和 分 级 。 

本 研究 利用 机 器 视觉 技术 对 柚子 轮廓 信息 进 
行 提取 ， 实 现 了 柚子 尺寸 、 果 形 等 参数 信息 的 无 
损 检测 ， 后 续 可 以 通过 以 下 途径 进一步 提高 柚子 
尺寸 检测 分 级 应 用 效果 。 

(1) 在 柚子 果 形 检测 分 类 判别 方面 ， 将 更 多 
深度 学 习 算 法 ， 如 支持 向 量 机 法 等 ， 应 用 到 柚子 
果 形 的 判别 ， 比 较 不 同 算法 的 优 劣 ， 寻 找 更 加 优 
越 的 判别 方法 ， 以 期 提高 柚子 果 形 的 识别 率 。 

(2) 本 研究 通过 主 视 、 左 镜像 和 右 镜 像 三 幅 
图 像 分 别 进 行 柚子 果 形 判别 ， 后 续 研 究 可 以 通过 
融合 三 幅 图 像 信 息 以 进行 柚子 果 形 的 判别 。 
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Detection and Grading Method of Pomelo Shape Based 
on Contour Coordinate Transformation and Fitting 
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Abstract: Automatic grading method of pomelo fruit according to the shape and size is urgently needed in the industry since the 
work mainly depends on artificial judgment currently. In this research, a method, which detected the vertical and horizontal size 
of pomelo by using contour coordinate transformation fitting, fruit shape feature extraction and direction angle compensation al- 
gorithm, while it determined the shape defects based on fruit shape index, was proposed. The image acquisition system was self- 
designed and built up with a CMOS camera, a dot matrix LED light source, a plane mirror, the computer, a box and brackets. 
The image data containing whole surface information of Shatian pomelo samples with different sizes and shapes were collected 
by this system. The G-B component grayscale image was chosen for denoising and segmentation. The Laplacian edge detection 
algorithm was implemented to extract the edge pixels of the fruit. The polynomial fitting method was applied to converse the 
rectangular coordinates to polar coordinates so that the fruit shape description was simplified. The characteristic point polar an- 
gle value was used to compensate the random direction of the vertical and horizontal diameters of the sample. Then the vertical 
and horizontal diameters of fruit were calculated after classifying the sample shapes into the spherical and the pear-like catego- 
ries. For the involved 168 pomelo samples, the average error, maximum absolute error and average relative error of the vertical 
diameters were 2.23 mm, 7.39 mm and 1.6% respectively, while these parameters of the horizontal diameters were 2.21 mm, 
7.66 mm and 1.4% respectively. The fruit shape discriminant model was established by using BP neural network algorithm 
based on the seven features extracted from the fitting function and verified by independent validation set including 3 peak 
heights, 3 peak widths and 1 trough value difference. The total recognition rate of shape identification was 83.7%. The results il- 
lustrated that the method had the potential to measuring the pomelo size and shape for grading fast and non-destructively. 

Key words: pomelo contour; fruit shape detection; back propagation neural network; coordinate system conversion; image pro- 


cessing; fruit shape discriminant model 
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